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Abstract 

We explore a connection of the deterministically forced Burgers equation 
for local velocity fields with probabilistic solutions of the Schrodinger bound- 
ary data problem. An issue of deducing the most likely interpolating dynamics 
from the given initial and terminal probability density data is here investi- 
gated to give account of the perturbation by external electromagnetic fields. 
A suitable extension of the Hopf-Cole logarithmic transformation allows to 
deal with nonconservative drift fields and implies the validity of appropriately 
generalised heat equations, which completely determine the dynamics. 
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1 Burgers velocity fields and related transport 
processes 

The Burgers equation with (usually without, [[[], |J) the forcing term F(x,t): 

d t v + (vV)v = uAv + F(x,t) (1) 

and especially its statistically relevant curlv = solutions, recently have acquired 
a considerable popularity in the variety of physical contexts. They range from an 
astrophysical issue of the stratified large-scale distribution of matter in the early 
Universe, f|, through accoustic turbulence dealing with intense noise in com- 
pressible liquids and gases, |J, to primitive fluid turbulence modeling in terms of 
the statistics of Burgers shocks in the low viscosity regime (enhanced by random 
initial data), [[3], [I], 0, §], eventually ending with the analysis of a fully developed 
Burgers turbulence that is regarded as a result of random forcing (stirring) of respec- 



tive velocity fields, |], [10], [TH| . It also pertains to the turbulence- without-pressure 



models, |12|], description of directed polymers in a random potential, fill, random 



interface growth problem, [13], and fluctuations/dispersion in deterministic or ran- 
dom flows, [[13], [14|, [l^, |l6j. Last but not least, all that implies quite nontrivial 
mathematics (very prohibitive with respect to the widely spread formal white noise 
manipulations), |17|, [18], |19| . 

Although Burgers velocity fields can be analysed on their own, frequently one 
needs a supplementary insight into the matter transport dynamics that is consistent 
with the chosen (Burgers) velocity field evolution. Then, the passive scalar (tracer 
or contaminant) advection-in-a-flow problem, ]14|, IT, [16| naturally appears through 
the parabolic dynamics: 

d t T + (£V)T = u/\T (2) 

Here, the Burgers velocity field defines the flow which amplifies the "passive scalar 
dispersion" (in analogy to the standard turbulent transport effects). This in turn 
implies the time-evolution of the concentration T(x,t) of the contaminant in the 
driving medium (incompressible fluid, Vv = 0, but infinitely compressible gas as 
well). While looking for the stochastic implementation of the microscopic (molecu- 



lar) dynamics, [^1], IT, [TJJ , it is assumed that the " diffusing scalar" (contaminant in 



the lore of early statistical turbulence models) obeys an Ito equation: 

dX{t) = v(x, t)dt + \f2vdW{t) (3) 
X(0) = x X(t) = x 
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where the given forced Burgers velocity field is perturbed by the noise term rep- 
resenting a molecular diffusion. In the Ito representation of diffusion-type random 

— * , — # 

variable X(t) one explicitly refers to the Wiener process \/2uW(t), instead of the 
usually adopted formal white noise integral f ff(s)ds, coming from the Langevin- 
type version of (3). 

Under these premises, we cannot view equations (l)-(3) as completely indepen- 
dent (disjoint) problems: the velocity field v cannot be arbitrarily inferred from (1) 
or any other velocity-defining equation without verifying the consistency conditions, 
which would allow to associate (2) and (3) with a well defined random dynamics 
(stochastic process), and Markovian diffusion in particular, p2 , [23| . In connection 



with the usage of Burgers velocity fields (with or without external forcing) which 
in (3) clearly are intended to replace the standard forward drift of the would-be- 
involved Markov diffusion process, we have not found in the literature any attempt 
to resolve apparent contradictions arising if (2) and/or (3) are defined by means of 

(!)• 

Moreover, rather obvious hesitation could have been observed in attempts to 
establish the most appropriate matter transport rule, if (l)-(3) are adopted. De- 
pending on the particular phenomeno logical departure point, one either adopts the 
standard continuity equation, || |j, that is certainly valid to a high degree of ac- 
curacy in the low viscosity limit v \ of (l)-(3), but incorrect on mathematical 
grounds if there is a diffusion involved and simultaneously a solution of (1) stands 
for the respective current velocity of the flow: d t p(x,t) = —V[v(x,t)p(x,t)} . Al- 
ternatively, following the white noise calculus tradition telling that the stochastic 
integral X(t) = f v(X(s), s)ds + f f](s)ds necessarily implies the Fokker-Planck 
equation, one adopts, d t p(x,t) = vAp(x,t) — V[v(x,t)p(x,t)] which is clearly 
problematic in view of the classic Mc Kean's discussion of the propagation of chaos 



for the Burgers equation, [2J, |25|, [26[ and the derivation of the stochastic " Burgers 



process" in this context: "the fun begins in trying to describe this Burgers motion 



as the path of a tagged molecule in an infinite bath of like molecules" , [^J . 

Moreover, an issue of the necessary correlation (cf. [TO, Chap. 7. 3, devoted to 
the turbulent transport and the related dispersion of contaminants) between the 
probabilistic Fokker-Planck dynamics of the diffusing tracer, and this of the passive 
tracer (contaminant) concentration (2), usually is left aside in the literature. 

To put things on the solid ground, let us consider a Markovian diffusion process, 
which is characterised by the transition probability density (generally inhomoge- 
neous in space and time law of random displacements) p(y, s, x, t) , < s < t < T, 
and the probability density p(x,t) of its random variable X(t) , < t < T. The 
process is completely determined by these data. For clarity of discussion, we do not 
impose any spatial boundary restrictions, nor fix any concrete limiting value of T 
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which, in principle, can be moved to infinity. 
The conditions valid for any e > 0: 

(a) there holds Um tls -±- s f\^ >e p(y, s,x,t)d 3 x = 0, 

(b) there exists a (forward) drift b(x,s) = lim t i s ^ f\£_g\ <e (y — x)p(x, s,y,t)d 3 y, 

(c) there exists a diffusion function (in our case it is simply a diffusion coefficient v) 
a(x,s) = Um tlsI ^ I\g-g\< e (y ~ x) 2 p(x, s,y,t)d 3 y, 

are conventionally interpreted to define a diffusion process, P3L E2I. Under suitable 



restrictions (boundedness of involved functions, their continuous differentiability) 
the function: 

g(x, s) = E{g(X(T))\X(s) = x; s < T} = J p(x, s, y, T)g(y, T)d 3 y (4) 

satisfies the backward diffusion equation (notice two sign changes in comparison 
with (2)) 

-d s g(x, s) = uAg(x, s) + [b(x, s)V\g(sc, s) . (5) 

Let us point out that the validity of (5) is known to be a necessary condition for the 
existence of a Markov diffusion process, whose probability density p(x, t) is to obey 
the Fokker-Planck equation (the forward drift b(x, t) replaces the previously utilized 
Burgers velocity v(x, t)): 

d t p(x, t) = uAp(x, t) - V[6(x, t)p(x , t)] (6) 

The case of particular interest, in the traditional nonequilibrium statistical physics 
literature, appears when p(y, s, x, t) is a fundamental solution of (5) with respect to 
variables y, s, \ 27[ g3[ , see however |f28| for an analysis of alternative situations. 



Then, the transition probability density satisfies also the second Kolmogorov (e.g. 
the Fokker-Planck) equation in the remaining x, t pair of variables. Let us empha- 
size that these two equations form an adjoint pair of partial differential equations, 
referring to the slightly counterintuitive for physicists, though transparent for math- 



ematicians, |3(], |3T], |32], |33], [33), [35], issue of time reversal of diffusions. 

After adjusting (3) to the present context, X(t) = J * b(X(s), s) ds + y/2uW(t) we 



can utilize the standard rules of the Ito stochastic calculus, [^, [33], [34], |55]| , to realise 
that for any smooth function f(x,t) of the random variable X(t) the conditional 
expectation value: 

lim Atl0 -^{Jp(x,t,y,t + At)f(y,t + At)d 3 y - f(x,t)] = (D + f)(X(t),t) = (7) 

= (d t + bV + uA)f(x,t) , 
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where X(t) = x, determines the forward drift b(x,t) (if we set components of X 
instead of /) and allows to introduce the local field of (forward) accelerations asso- 
ciated with the diffusion process, which we constrain by demanding (see e.g. Refs. 



36| , p4| , |35| for prototypes of such dynamical constraints): 



(D 2 + X)(t) = (D + b)(X(t),t) = (d t b+(bV)b + vAb)(x,t) = F(x,t) (8) 

where X(t) = x and, at the moment arbitrary, function F(x,t) may be interpreted 
as an external forcing applied to the diffusing system, ||29|| . 

In particular, if we assume that drifts remain gradient fields, curl 6 = 0, under the 
forcing, then those that are allowed by the prescribed choice of F(x, t) must fulfill the 
compatibility condition (notice the conspicuous absence of the standard Newtonian 



minus sign in this analogue of the second Newton law, see e.g. |37| , |34| , |29| ) 



p(x, t) = Vfi(x , t) (9) 
which establishes the Girsanov-type martingale connection of the forward drift 



b(x,t) = 2z/V$(x, t) with the Feynman-Kac, cf. [29, 28], potential Q(x,t) of the 
chosen external force field. 

One of the distinctive features of Markovian diffusion processes with the pos- 
itive density p(x, t) is that the notion of the backward transition probability den- 
sity p*(y,s,x,t) can be consistently introduced on each finite time interval, say 
< s < t < T: 

p(x,t)p*(y,s,x,t) = p(y,s,x,t)p(y,s) (10) 

so that / p(y, s)p(y, s, x,t)d 3 y = p(x,t) and p(y,s) = J p*(y, s,x,t)p(x,t)d 3 x. This 
allows to define the backward derivative of the process in the conditional mean (cf. 
29], |38"| [39| for a discussion of these concepts in case of the most traditional Brownian 



motion and Smoluchowski-type diffusion processes) 

Hm Atl0 -^[x- j p*(y,t- At,x,t)yd 3 y] = (D-X)(t) = t(X(t), t) (11) 

(D_f)(X(t),t) = (d t + lV-uA)f(X(t),t) 
Accordingly, the backward version of the acceleration field reads 

(D*_X)(t) = (DlX)(t) = F(X(t),t) (12) 
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where in view of 6* = b — 2u'Vlnp we have explicitly fulfilled the forced Burgers 
equation (cf. (1)): 

OA + (6*V)6* - 1/A6, = F (13) 

and, [[33], |34], |29j, under the gradient-drift field assumption, curl = 0, we deal with 
F(x, t) given by (9) again. A notable consequence of the involved backward Ito 
calculus is that the Fokker-Planck equation (6) can be transformed to an equivalent 
form of: 

d t p(x,t) = -uAp(x,t) - V[t(x,t)p(x,t)] (14) 

with the very same initial (Cauchy) data po(x) = p(x,0) as in (6). 

At this point let us recall that equations (5) and (6) form a natural adjoint pair 
of equations for the diffusion process in the time interval [0, T\. Clearly, an adjoint 
of (14) reads: 

d s f(x, s) = uAf(x, s) - [l(x, s)V]/(f , s) (15) 

where: 

f(x,s) = / p*(y,Q,x,s)f(y,Q)d 3 y , (16) 



to be compared with (4), (5) and the previously mentioned passive scalar dynamics 
(2), see e.g. also |21]. Here, manifestly, the time evolution of the backward drift is 
governed by the Burgers equation, and the diffusion equation (15) is correlated (via 
the definition (10)) with the probability density evolution rule (14). 

This pair only can be consistently utilized if the diffusion proces is to be driven 
by forced (or unforced) Burgers velocity fields. 

On the other hand, the study of diffusion driven by the Burgers flow may begin 
from first solving the Burgers equation (12) for a chosen external force field, next 
specifying the probability density (14), eventually ending with the corresponding 
"passive contaminant" concentration dynamics (15), (16). All that remains in a 
perfect agreement with the heuristic discussion of the concentration dynamics given 
in Ref. [IE|, Chap. 7.3. where the "backward dispersion" problem with "time 



running backwards" was found necessary to predict the concentration. 

Let us notice that the familiar logarithmic Hopf-Cole transformation, [[|, |43| , 
of the Burgers equation into the generalised diffusion equation (yielding explicit 
solutions in the unforced case) has received a generalisation in the framework of the 
so called Schrodinger interpolation problem, p|, JJH 123, Ejl, fM fUJfl, see also Wl 



In its recent reformulation, , the solution in terms of the interpolating Markovian 



diffusion process is found to rely on the adjoint pairs of parabolic equations, like e.g. 
(5), (6) or (14), (15). In case of gradient drift fields the process can be determined 
by checking (this imposes limitations on the admissible force field potential) whether 
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the Feynman-Kac kernel 

k(y, s, x,t) = J exp[- J c(w(r), T)d,T]dfj,^(u) (17) 
is positive and continuous in the open space-time area of interest (then, additional 



limitations on the path measure need to be introduced, p9|, (£!]), and whether it 



gives rise to positive solutions of the adjoint pair of generalised heat equations: 

d t u(x, t) = vAu(x, t) — c(x, t)u(x, t) (18) 

d t v(x, t) = —vAv(x, t) + c(x, t)v(x, t) 

Here, a function c(x,t) = j^Q(x,t) is restricted only by the positivity and conti- 
nuity demand for the kernel (17), see e.g. fZS). In the above, d^% s l{uj) is the 
conditional Wiener measure over sample paths of the standard Brownian motion. 

Solutions of (18), upon suitable normalisation give rise to the Markovian diffusion 
process with the factorised probability density p(x,t) = u(x,t)v(x,t) wich interpo- 
lates between the boundary density data p(x, 0) and p(x,T), with the forward and 
backward drifts of the process defined as follows: 

r, ^ \ Vv(x, t) , . 

^•^W (19) 

u(x, t) 

in the prescribed time interval [0, T\. These formulas imply that the compatibil- 
ity condition (9), connecting the drifts with the a priori chosen function c(x,t) = 
-^Vt(x,t), can be easily established. 

The transition probability density of this process reads: 

p(y, s, x, t) = k(y, s, x, t) V ^ S . (20) 

while the corresponding (since p{x,t) is given) transition probability density, (10), 
of the backward process has the form: 

uiv S) 

P*(y,s,x,t) = k(y,s,x,t) ' . (21) 

u{x, t) 



Obviously, f28" , [3^| , in the time interval < s < t < T there holds: 

u{x , t)= I u {y)k(y, s, x, t)d 3 y (22) 



v(y,s) = / k(y,s,x,T)v T (x)d 3 x . 



By defining = logu, we immediately recover the traditional form of the 
Hopf-Cole transformation for Burgers velocity fields: = — 2z/V < £*. In the special 
case of the standard free Brownian motion, there holds b(x,t) = while b*(x,t) = 
—2uVlog p(x, t). 

2 The problem of electromagnetic forcing in the 
Schrodinger interpolating dynamics 

It turns out the the crucial point of our previous discussion lies in a proper choice of 
the strictly positive and continuous, in an open space-time area, function k(y, s, x, t) 
which, if we wish to construct a Markov process, has to satisfy the Chapman- 
Kolmogorov (semigroup composition) equation. It has led us to consider a pair 
of adjoint partial differential equations, (18), as an alternative to either (5), (6) or 
(14), (15). 

In the quantally oriented literature dealing with Schrodinger operators and their 
spectral properties, ||6j [47|, pj] , the potential c(x, t) in (17), (18) is usually assumed 



to be a continuous and bounded from below function, but these restrictions can be 
substantially relaxed (unbounded functions are allowed in principle) if we wish to 
consider general Markovian diffusion processes and disregard an issue of the bound 
state spectrum and this of the ground state of the (selfadjoint) semigroup generator, 



22 1 . Actually, what we need is merely that properties of c(x, t) allow for the 



kernel k, (17), which is positive and continuous. Taking for granted that suitable 



conditions are fulfilled, |28| , we can immediately associate with equations (18) 
an integral kernel of the time-dependent semigroup (the exponential operator should 
be understood as time-ordered expression, since in general H(t) may not commute 
with H(t') for t^t'): 

k(y,s,x,t) = [exp(- f H(r)dr)](y,x) (23) 

J s 

where H(t) = —uA + c(r) is the pertinent semigroup generator. Then, by the 
Feynman-Kac formula, E2J, we get an expression (17) for the kernel, which in turn 
yields (19)-(22), see e.g. p8| . As mentioned before, (20) combined with (17) sets a 
Girsanov probabilistic (martingale) connection between the Wiener measure corre- 
sponding to the standard Brownian motion with b(x, t) = 0, and this appropriate 
for the diffusion process with a nonvanishing drift b(x, t), curl 6 = 0. 



The above formalism is known, [29], to encompass the standard Smoluchowski- 



type diffusions in conservative force fields. 
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Strikingly, an investigation of electromagnetically forced diffusions has not been 
much pursued in the literature, although an issue of deriving the Smoluchowski- 
Kramers equation (and possibly its large friction limit) from the Langevin-type 
equation for the charged Brownian particle in the general electromagnetic field has 
been relegated in Ref. [fUE] , Chap. 6.1 to the status of the innocent-looking exercise 



(sic !). On the other hand, the diffusion of realistic charges in dilute ionic solutions 
creates a number of additional difficulties due to the apparent Hall mobility in terms 
of mean currents induced by the electric field (once assumed to act upon the system), 
see e.g. @ [50| and 



In connection with the electromagnetic forcing of diffusing charges, the gradient 
field assumption imposes a severe limitation if we account for typical (nonzero circu- 
lation) features of the classical motion due to the Lorentz force, with or without the 
random perturbation component. The purely electric forcing is simpler to handle, 
since it has a definite gradient field realisation, see e.g. for a recent discussion 
of related issues. The major obstacle with respect to our previous (Section 1) dis- 
cussion is that, if we wish to regard either the force F, (8), (12), or drifts b, 6* to 
have an electromagnetic origin, then necessarily we need to pass from conservative 
to non-conservative fields. This subject matter has not been significantly exploited 
so far in the nonequilibrium statistical physics literature. 

Interestingly, in the framework of the Onsager-Machlup approach towards an 
identification of most probable paths which one can associate with the underlying 
diffusion process, 0, [54], [55], the non-conservative model system has been inves- 
tigated, |p6[ . Namely, an effectively two-dimensional Brownian motion was ana- 
lyzed, whose three-dimensional forward drift b(x), 6 3 = in view of d x b\ ^ d y b 2 
has curlb ^ 0. Then, by the standard variational argument with respect to the 



Wiener- Onsager-Machlup action, [54, Bq 



1 ft 2 1 

I{L&x,t);t 1 ,t 2 } = — {-[S-b(x,t)} 2 + uVb(x,t)}dt , (24) 

LV Jt! I 

the most probable trajectory, about which major contributions from (weighted) 
Brownian paths are concentrated, was found to be a solution of the Euler-Lagrange 
equations, which are formally identical to the equations of motion 

'q cl = E + <j d xB (25) 

of a classical particle of unit mass and unit charge moving in an electric field E and 
the magnetic field B. The electric field (to be compared with Eq. (9)) is given by: 

E = -V$ (26) 
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while the magnetic field has the only nonvanishing component in the z-direction of 

B = curl b = {0, 0, d x b 2 - dyh} , (27) 

— * — * — * — * 

Clearly, B = curl A where A=b is the electromagnetic vector potential. The simplest 
example is a notorious constant magnetic field defined by b\{x) = — ^x 2 , b 2 (x) = 
fa*. 

One immediately realizes that the Fokker-Planck equation in this case is incom- 
patible with traditional intuitions underlying the Smoluchowski-drift identification: 
the forward drift is not proportional to an external force, but to an electromag- 
netic potential. Nevertheless, the variational information drawn from the Onsager- 
Machlup Lagrangian involves the Lorentz force-driven trajectory. Hence, some prin- 
cipal effects of the electromagnetic forcing are present in the diffusing system, whose 
drifts display an "unphysical" (gauge dependent) form. 

On the other hand, if we accept this "unphysical" random motion to yield the 
representation dX(t) = A(X(t),t)dt + y/2udW(t), supplemented by the correspond- 
ing pair (5), (6) of adjoint diffusion equations with A(x,t) replacing b(x,t), then 
standard rules of the Ito stochastic calculus, (8), tell us that 

d2 

(D 2 + X)(t) =d t A+ (AV)A + vAA = -— {x 1: x 2 ,0} = -E(x) (28) 

where E(x) = ^-{xi, x 2 ,0}, if calculated from (26). 

We thus arrive at the purely electric forcing with reversed sign (if compared with 
that coming from the Onsager-Machlup argument, (26)) and somewhat surprisingly, 
there is no impact of the previously discussed magnetic motion on the level of dy- 
namical constraints (8), (13). The adopted recipe is thus incapable of producing the 
magnetically forced diffusion process that conforms with arguments of Section 1. 

In below, we shall follow more abstract route, based on invoking the Feynman- 
Kac kernel idea (23), f28fl . This approach has a clear advantage of elucidating the 
generic issues that hamper attempts to describe electromagnetically forced diffusion 
processes. 

Usually, the selfadjoint semigroup generators attract the attention of physicists 
in connection with the Feynman-Kac formula. A typical route towards incorporating 
electromagnetism comes from quantal motivations via the minimal electromagnetic 
coupling recipe which preserves the selfadjointness of the generator (Hamiltonian 
of the system). As such, it constitutes a part of the general theory of Schrodinger 
operators. A rigorous study of operators of the form — A + V has become a well 
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developed mathematical discipline, |4"6|| . The study of Schrodinger operators with 
magnetic fields, typically of the form —(V — iA) 2 + V , is less advanced, although 
specialised chapters on the magnetic field issue can be found in monographs devoted 



to functional integration methods, [46, 57], mostly in reference to seminal papers 
§0. 

From the mathematical point of view, it is desirable to deal with magnetic fields 
that go to zero at infinity, which is certainly acceptable on physical grounds as well. 
The constant magnetic field (see e.g. our previous considerations) does not meet 
this requirement, and its notorious usage in the literature makes us (at the moment) 
to decline the asymptotic assumption and inevitably fall into a number of serious 
complications. 

One obvious obstacle can be seen immediately by taking advantage of the ex- 
isting results, [5Ej|. Namely, an explicit expression for the Feynman-Kac kernel in a 



constant magnetic field, introduced through the the minimal electromagnetic cou- 
pling assumption H(A) = — |(V — iA) 2 , is available (up to irrelevant dimensional 



constants): 



eipW -) 1(g ,a = [ i — g 1( JLr (29) 



1 B B B 

xexp{-—(x 3 - y 3 ) 2 - —coth(—t)[(x 2 - y 2 f + (xi - yi) 2 } - i—{x x y 2 - x 2 yi)} . 

Clearly, it is not real (hence non-positive and directly at variance with the major 
demand in the Schrodinger interpolation problem, as outlined in Section 1), except 
for directions y that are parallel to a chosen x. 

Consequently, a bulk of the well developed mathematical theory is of no use for 
our purposes and new techniques must be developed for a consistent description of 
the electromagnetically forced diffusion processes along the lines of Section 1, i.e. 
within the framework of Schrodinger's interpolation problem. 



3 Forcing via Feynman-Kac semigroups 

The conditional Wiener measure d^% s 2{oj), appearing in the Feynman-Kac kernel 
definition (17), if unweighted (set c(uj(t),t) = 0) gives rise to the familiar heat 
kernel. This, in turn, induces the Wiener measure Pw of the set of all sample 
paths, which originate from y at time s and terminate (can be located) in the Borel 
set A e -R 3 after time t — s: P\y[A} = J A d 3 x J d/jr^hu) = J A d/j where, for simplicity 

of notation, the (y, t — s) labels are omitted and stands for the heat kernel. 
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Having denned an Ito diffusion X(t) = Jq b(x, u)du + y/2uW(t) we are interested 
in the analogous path measure: -P^[^4] = J A dx J diJ& s }(ujg) = J A du(X). 

Under suitable (stochastic, |^J) integrability conditions imposed on the forward 
drift, we have granted the absolute continuity Px <C Pw of measures, which implies 
the existence of a strictly positive Radon-Nikodym density. Its canonical Cameron- 
Martin-Girsanov form, EoT, E3, reads: 



d ^ X \y,s,x,t) = exp^-[f t b(X{u),u)dX(u) - \ f \b{X{u),u)] 2 du} . (30) 



dfx 2v Js 2 

If we assume that drifts are gradient fields, curl 6 = 0, then the Ito formula 
allows to reduce, otherwise troublesome, stochastic integration in the exponent of 
(30), |46], [57]], to ordinary Lebesgue integrals: 



fb(X{u)iu)dX{u) = <Z>(X(t),t) - [ t du[d t $ + l-Vb](X(u),u) . 
2v Js Js 2 

(31) 

After inserting (31) to (30) and next integrating with respect to the conditional 
Wiener measure, on account of (9) we arrive at the standard form of the Feynman- 
Kac kernel (17). Notice that (31) establishes a probabilistic basis for logarithmic 
transformations (19) of forward and backward drifts: b = 2u'Vlogv = 2z/V$, 6* = 
—2vVlog u = —2^X7$*. The forward version is commonly used in connection with 
the transformation of the Fokker-Planck equation into the generalised heat equation, 



37]. 0, |29|. The backward version is the Hopf-Cole transformation, mentioned in 
Section 1, used to map the Burgers equation into the very same generalised heat 
equation as in the previous case, 0, f45f . 

— * 

However, presently we are interested in non-conservative drift fields, curl b ^ 0, 
and in that case the stochastic integral in (30) is the major source of computational 



difficulties, HB], I57L B3|, for nontrivial vector potential field configurations. It explains 



the virtual absence of magnetically forced diffusion problems in the nonequilibrium 
statistical physics literature. 

At this point, some steps of the analysis performed in Ref. ]63|] in the context 



of the "Euclidean quantum mechanics", cf. also p5| , are extremely useful. Let us 
emphasize that electromagnetic fields we utilize, are always meant to be ordinary 
Maxwell fields with no Euclidean connotations (see e.g. Chap. 9 of Ref. |57j for the 
Euclidean version of Maxwell theory). 

Let us consider a gradient drift-field diffusion problem according to Section 1, 
with (17), (31) involved and thus an adjoint pair (18) of parabolic equations com- 
pletely defining the Markovian diffusion process. Furthermore, let A(x) be the 
time-independent vector potential for the Maxwellian magnetic field B = curl A. 
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We pass from the gradient realisation of drifts to the new one, generalizing (19), 
for which the following decomposition into the gradient and nonconservative part is 
valid: 

b{x , t) = 2v log t) - A(x) , (32) 

We denote 9(x,t)=exp[$>(x,t)] and admit that (32) is a forward drift of an Ito 
diffusion process with a stochastic differential dX(t) = — A]dt + \/2udW{t). 

On purely formal grounds, we deal here with an example of the Cameron-Martin- 
Girsanov transformation of the forward drift of a given Markovian diffusion process 
and we are entitled to ask for a corresponding measure transformation, (30). 

To this end, let us furthermore assume that 9(x, t) = 9 solves a partial differential 
equation 

d t 9 = -v\y - ^-A(x)} 2 9 + c{x, t)9 (33) 

with the notation c(x,t) = ^Q(x,t) patterned after (9). Then, by using the Ito 
calculus and (32), (33) on the way, see e.g. Ref. p3]|, we can rewrite (30) as follows: 



^±(y,s,x,t)= (34) 

= ex V^\\\^ ~ A)(X(u),u)dX(u) - l - J\2v™ - A) 2 (X(u),u) du] 
6 ^ X ^' t \exp[~— f {A(u)dX(u) + is(VA)(X(u))du + n(X(u),u)du}} , 



9(X(s),s) 2u 

where X(s) = y,X(t) = x. 

More significant observation is that the Radon-Nikodym density (34), if inte- 
grated with respect to the conditional Wiener measure, gives rise to the Feynman- 
Kac kernel (23) of the non- self adjoint semigroup (suitable integrability conditions 
need to be respected here as well, p3[|), with the generator = — z/[V— -^A(x)} 2 + 
c(x, t) defined by the right-hand-side of (33): 

d t 6(x,t) = H A 9(x,t) = 

= [-i/A + A(x)V + l -(yA(x)) - ^[A(x)f + c(x,t)]9(x,t) (35) 
= -vA9(x, t) + A(x)V9(x, t)+c A (x, t)9(x, i) . 

Here: 

c A (x,t) = c(x,t) + ~(VA)(x) - ^\A(x)} 2 . (36) 
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An adjoint parabolic partner of (35) reads: 



d t 9, = -H* A 9, = vA6« + V[A{x)6*\ - c A (x, t)0* = (37) 

= z/[V + ^A{x)] 2 6* - c(x, t)6, . 

Consequently, our assumptions (32), (33) involve a generalization of the adjoint 
parabolic system (18) to a new adjoint one comprising (33), (37). Obviously, the 
original form of (18) is immediately restored by setting A = 0, and executing obvious 
replacements 6* — > u, 6 — > v. 

Let us emphasize again, that in contrast to Ref. where the non-Hermitean 



generator 2vH^, (33), has been introduced as "the Euclidean version of the Hamil- 
tonian" H = — 2u 2 (V — ^A) 2 + Q, our electromagnetic fields stand for solutions of 
the usual Maxwell equations and are not Euclidean at all. 

As long as the coefficient functions (both additive and multiplicative) of the 
adjoint parabolic system (35), (37) are not specified, we remain within a general 
theory of positive solutions for parabolic equations with unbounded coefficients (of 
particular importance, if we do not impose any asymptotic fall off restrictions), 
]2"7| , ISO , |HT|, p3|. The fundamental solutions, if their existence can be granted, usually 
live on space-time strips, and generally do not admit unbounded time intervals. We 
shall disregard these issues at the moment, and assume the existence of fundamental 
solutions without any reservations. 

By exploiting the rules of functional (Malliavin, variational) calculus, under an 
assumption that we deal with a diffusion (in fact, Bernstein) process associated 
with an adjoint pair (35), (36), it has been shown in Ref. |63[] that if the forward 
conditonal derivatives of the process exist, then (D+X)(t) = - A — b(x,t), 

(32) and: 

{D\X)(t) = {D + X)(t) x curlA(x) + Vfi(f,t) + vcurl (curl A(x)) (38) 

where X(0) = 0, X{t) = x, x denotes the vector product in R 3 and 2uc = Q. 

Since B = curl A = fioH, we identify in the above the standard Maxwell equation 
for curl H comprising magnetic effects of electric currents in the system: curl B = 

fio{D + o~oE + J ext ] where D = e^E while J ext represents external electric currents. In 
case of E = 0, the external currents only would be relevant. A demand curl curl A = 
V(VA) — A A = corresponds to a total absence of such currents, and the Coulomb 
gauge choice VA = would leave us with harmonic functions A(x). 

Consequently, a correct expression for the magnetically implemented Lorentz 
force has appeared on the right-hand-side of the forward acceleration formula (38), 
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with the forward drift (32) replacing the classical particle velocity q of the classical 
formula (25). 

The above discussion implicitly involves quite sophisticated mathematics, hence 
it is instructive to see that we can bypass the apparent complications by directly 
invoking the universal definitions (7) and (11) of conditional expectation values, 
that are based on exploitation of the Ito formula only. Obviously, under an assump- 
tion that the Markovian diffusion process with well defined transition probability 
densities p(y, s, x, t) and p*(y, s, x, t), does exist. 

We shall utilize an obvious generalization of canonical definitions (19) of both 
forward and backward drifts of the diffusion process defined by the adjoint parabolic 
pair (18), as suggested by (32) with A = A(x): 

b = 2I/-V - A , t = - A . (39) 

We also demand that the corresponding adjoint equations (35), (37) are solved by 
8 and 6* respectively. 

Taking for granted that identities (D + X)(t) = b(x, t), X(t) — x and (D-X)(t) = 
b*(x,t) hold true, we can easily evaluate the forward and backward accelerations 
(substitute (39), and exploit the equations (35), (37)): 

(D+b)(X(t),t) = d t b+ (W)b + uAb= (40) 

= b x B + v curl B + V VL 

and 

(DX)(X(t),t) = d t t + (&*V)6* - uAK = (41) 
= K x B - v curl B + V VL . 

Let us notice that the forward and backward acceleration formulas do not coin- 
cide as was the case before ( cf. Eq. (8), (12)). There is a definite time-asymmetry 
in the local description of the diffusion process in the presence of general magnetic 
fields, unless curl B = 0. The quantity which is explicitly time-reversal invariant 
can be easily introduced: 

v(x,t) = -(b + &)(x,t)=> (42) 

±(Dl + D 2 _)(X(t))=vxB* + Vn . 

As yet there is no trace of Lorentzian electric forces, unless extracted from the term 
VQ(x,t). This step we shall accomplish in Section 4. 
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For a probability density 6*6 = p of the related Markovian diffusion process, 
34, 28 1, we would have fulfilled both the Fokker-Planck and the continuity equations: 



d t p = z/Ap — V(bp) = —V(vp) = —uAp — V(6*p), as before (cf. Section 1). 

In the above, the equation (41) can be regarded as the Burgers equation with 
a general external magnetic (plus other external force contributions if necessary) 
forcing, and its definition is an outcome of the underlying mathematical structure 
related to the adjoint pair (33), (37) of parabolic equations. Our construction shows 
that the solution of the magnetically forced Burgers equation needs to be sought in 
the form (39). 



4 Schrodinger's interpolation in a constant mag- 
netic field and quantally inspired generalisations 

Presently, we shall confine our attention to the simplest case of a constant magnetic 
field, defined by the vector potential A = { — —#2, +§ %i, 0}. Here, B = {0,0, B], 
VA = and curl B = which significantly simplifies formulas (32)-(42). 

As emphasized before, most of our discussion was based on the existence assump- 
tion for fundamental solutions of the (adjoint) parabolic equations (33), (37). For 
magnetic fields, which do not vanish at spatial infinities (hence for our "simplest" 
choice), the situation becomes rather complicated. Namely, an expression for 

B 2 

cj(f , t) = c(f, t) - —{x\ + x\) (43) 
Ibis 

includes a repulsive harmonic oscillator contribution. 

For the existence of a well defined Markovian diffusion process it appears neces- 
sary that a nonvanishing contribution from an unbounded from above c(x, t) would 
counterbalance the harmonic repulsion. To see that this must be the case, let us for- 
mally constrain 6(x,t) = exp(<&(x,t)) to yield (in accordance with (9)) the identity: 

c (x, t) = dt$ + z/[V$] 2 + z/A$ = (44) 

Then, we deal with the simplest version of the adjoint system (35), (37) where, in 
view of VA = = c, there holds: 

d t 6 = - ^A 2& = -" Ad + AV6 - ^[Af 6 (45) 
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With our choice, curl A = {0, 0, B}, the equations (45) do not possess a funda- 
mental solution, which would be well defined for all (x,t) G R 3 x R + : everything 
because of the harmonic repulsion term in the forward parabolic equation. In the 
Appendix, we shall prove that the function k(y, s,x,t), : 

k(y,s,x,t) = [- tttt? d(— r) 1/2 (46) 

IBB B 

x exp { ~ 2(t - s) ^ 3 ~ y ^ 2 ~~ ~l COt ^ ~2 ^ ~~ ^ ^ ^°° 2 ~ V2 ^ + ~ Vl ^ ^ ~ ~2 ( X lV2 ~~ X2Vl ^ ^ 

only when restricted to times t — s < is an acceptable example of a unique positive 
(actually, positivity extends to times t — s < ^) fundamental solution of the system 
(44), (rescaled to yield v — > |). Here, formally, (46) can be obtained from the 

— * —* 

expression (29) by the replacement A — > —iA. 

An immediate insight into a harmonic repulsion obstacle can be achieved after 
an x — y-plane rotation of Cartesian coordinates: x[ = XiCos(ut) — x 2 sin(ut), x' 2 = 
Xisin(ut) + X2cos(cot), x' 3 = x 3 , t' — t with cu = Then, equations (45) get 

transformed into an adjoint pair: 

fyO = -uA'6 - u 2 (x' 2 + x'i)9 (47) 

dt>6* = i/A'0, + ou 2 (x' 2 + xf)6* . 

Notice that the transformation oj — > iuj would replace repulsion in (47) by the 
harmonic attraction. On the other hand, we can get rid of the repulsive term by 
assuming that c(x, t), (43) does not identically vanish. For example, we can formally 
demand that, instead of (44), c(x,t) = +^;{x\ + x 2 ,) plays the role of an electric 
potential. Then, harmonic attraction replaces repulsion in the final form of equations 

(35), (37). ' 

As a byproduct, we are given a transition probability density of the diffusion 
process governed by the adjoint system (cf. (28)): 

d t 6 = -vA6 + AV9 (48) 

d t e* = vAe* + Ave, . 

with A = § { —X2,xi,0}. Namely, by means of the previous x — y-plane rotation, 
(48) is transformed into a pair of time adjoint heat equations: 

d t ,e = -vA'e , d t ,e* = v a% (49) 

whose fundamental solution is the standard heat kernel. 
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Finding explicit analytic solutions of rather involved equations (35), (37) is a 
formidable task on its own, in contrast to much simpler-unforced or conservatively 
forced dynamics issue. 

Interestingly, we can produce a number of examples by invoking the quantum 
Schrodinger dynamics. This quantum inspiration has been proved to be very useful 
in the past, |34|, j35fl . At this point, we shall follow the idea of Ref. p8| where 



the strategy developed for solving the Schrodinger boundary data problem has been 
applied to quantally induced stochastic processes (e.g. Nelson's diffusions, [[33, 



They were considered as a particular case of the general theory appropriate for 
nonequilibrium statistical physics processes as governed by the adjoint pair (18), 
and exclusively in conjunction with Born's statistical postulate in quantum theory. 

The Schrodinger picture quantum evolution is then consistently representable 
as a Markovian diffusion process. All that follows from the previously outlined 
Feynman-Kac kernel route, P3L EB, 133, I3DI, HDJ, based on exploiting the 



adjoint pairs of parabolic equations. However, the respective semigroup theory has 
been developed for pure gradient drift fields, hence without reference to any impact 
of electromagnetism on the pertinent diffusion process. And electromagnetism is 
definitely ubiquitous in the world of quantum phenomena. 

Let us start from an ordinary Schrodinger equation for a charged particle in an 
arbitrary external electromagnetic field, in its standard dimensional form. To con- 
form with the previous notation let us absorb the charge e and mass m parameters 
in the definition of A(x) and the potential 4>(x), e.g. we consider B instead of 
and instead <p/m. Additionally, we set v instead of Then, we have: 

id t if;(x, t) = -z/(V - ^Afifj{x, t) + —<f)(x)if)(x, t) . (50) 

The standard Madelung substitution if> = exp(R + iS) allows to introduce the 
real functions 9 = exp(R + S) and 9* = exp(R — S) instead of complex ones if), if). 
They are solutions of an adjoint parabolic system (35), (37), where the impact of 
(50) is encoded in a specific functional form of, otherwise arbitrary potential c(x, t): 

c(x,t) = ^Q(x,t) = ^[2Q(x,t) -0(f)] (51) 

Q(x,t) = 2z/ 2 A ^jf; t) = 2u 2 [AR(x,t) + [VR(x,t)] 2 ] . 

The quantum probability density p(x, t) = if>(x, t)if>(x, t) = 9(x, t)9*(x, t) displays 
a factorisation p = 99* in terms of solutions of adjoint parabolic equations, which we 
recognize to be characteristic for probabilistic solutions (Markov diffusion processes) 
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of the Schrodinger boundary data problem (cf. Section 1), [[34], |29], [28], [39]. It is easy 
to verify the validity of the Fokker-Planck equation whose forward drift has the form 
(38). Also, equations (40), (41) do follow with Q = 2Q - 0. 

By defining E = — V0 (with <p utilised instead of —<f>), we immediately arrive 
at the complete Lorentz force contribution in all acceleration formulas (before, we 
have used curl B = 0): 

d t b+ (6V)6 + z/A6 = 6 x B + E + ucurl B + 2VQ (52) 

d t h + (6* V)6* - vAh =K x B + E - ucurl B + 2VQ 

Moreover, the velocity field named the current velocity of the flow, v = 1(6 + 6*), 
enters the familiar local conservation laws (see also P9l for a discussion of how the 



"quantum potential" Q affects such laws in case of the standard Brownian motion 
and Smoluchowski-type diffusion processes) 

dp = -V(vp) (53) 

d t v + (vV)v = vx B + E + VQ . 

A comparison with (33)-(43) shows that equations (50)-(53) can be regarded 
as the specialized version of the general external forcing problem with an explicit 
electromagnetic (Lorentz force-inducing) contribution and an arbitrary term of non- 
electromagnetic origin, which we denote by c(x, t) again. Obviously, c is represented 
in(51),by±Q(x,t). 

We have therefore arrived at the following ultimate generalization of the adjoint 
parabolic system (18), that encompasses the nonequilibrium statistical physics and 
essentially quantum evolutions on an equal footing (with no clear-cut discrimination 
between these options, as in Ref. ||28|| ) and gives rise to an external (Lorentz) 
electromagnetic forcing: 

d t 6(x, t) = [-u(V - ^Af - + c(x, t)]6(x, t) (54) 

dMx,t) = [z/(V + ^A) 2 + -*-</>(£) - c(x,t)%(x,t) . 

A subsequent generalisation encompassing time-dependent electromagnetic fields is 
immediate. 

The adjoint parabolic pair (54) of equations can thus be regarded to determine a 
Markovian diffusion process in exactly the same way as (18) did. If only a suitable 
choice of vector and scalar potentials in (54) guarrantees a continuity and positivity 
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of the involved semigroup kernel (take the Radon-Nikodym density of the form (34), 
with Q — > —0 + Q , and integrate with respect to the conditional Wiener measure), 
then the mere knowledge of such integral kernel suffices for the implementation of 
steps (18)-(22), with u — > 9*, v — > 9. To this end it is not at all necessary that 
k(x,s,y,t) is a fundamental solution of (54). A sufficient condition is that the 
semigroup kernel is a continuous (and positive) function. The kernel may not even 
be different iable, see e.g. Ref. for a discussion of that issue which is typical for 
quantal situations. 

After adopting (54) as the principal dynamical ingredient of the electromagnet- 
ically forced Schrodinger interpolation, we must slightly adjust the emerging accel- 
eration formulas. Namely, they have the form (52) , but we need to replace 2Q(x, t) 
by, from now on arbitrary, potential Q(x,t) = 2vc{x,t). The second equation in 
(53) also takes a new form: 

d t v+ (vV)v = v x B + E + V(fl-Q) (55) 

see e.g. Ref. |29| for more detailed explanation of this step. The presence in (54) of 
the density-dependent — VQ term finds its origin in the identity b — b*= 2uVp(x, t) 
and is a necessary consequence of the involved (forced in the present case) Brownian 
motion, see e.g. [p4T BH]. 



Finally, the second equation (52) with f2 replacing 2Q is the most general form of 
the Burgers equation with an external forcing, where the electromagnetic (Lorentz 
force) contribution has been extracted for convenience. Solutions of this equa- 
tion must be sought for in the form (39), which generalizes the logarithmic Hopf- 
Cole transformation to non-gradient drift fields. Equations (54) are the associated 
parabolic partial differential (generalised heat) equations, that completely determine 
probabilistic solutions (Markovian diffusion processes) of the Schrodinger boundary 
data (interpolation) problem, for which the forced Burgers velocity fields play the 
role of backward drifts of the process. 



Appendix 

General criterions, allowing to decide when the semigroup kernel associated with the 
adjoint pair (54) is a continuous and positive function, belong to a specialised branch 
of mathematical theory of parabolic partial differential equations and mathematical 
physics. A list of restrictions suitable for our case can be found in Section 4 of Ref. 

0- 

To give a flavour of a typical argumentation, we shall present a more selective 
discussion of conditions under which the quantally inspired parabolic problem (33), 



20 



(37), (51) admits fundamental solutions. Let us point out that quite generally 
fundamental solutions do not exist (the kernels do !) in such case, cf. |28], ^0 . 
Let us rescale the problem (50) so that we are interested in: 

id t if){x, t) = --(V - iA) 2 ip{x, t) + (f)(x)if)(x, t) . (A.l) 

with the initial data if)o(x) G L 2 (R 3 ). We assume that A iy i = 1,2,3, G C l (R 3 ) i.e. 
have continuous first derivatives, while <p(x) G C(R 3 ) and is bounded from below, 
e.g. 4>(x) > —M for some positive M. It is well known that under these restrictions 
the operator H = |(iV + A) 2 + <p is essentially selfadjoint on C^°(R 3 ), i.e. the space 
of smooth functions with a compact support, Thus ipt = exp(—itH)ip is a 

solution of (A.l) for every time t and any initial condition ipo G D(H). Here, if 
denotes the closure of an operator H. 

Let us fix some T > and let us assume that if) (x, t) ^ for all x E R 3 and 
t G [0, T]. Then, the Madelung substitution if) = exp(R + iS) implies: 

d t R = -VRVS - ^AS + -VA + AVR (A2) 

d t S = ^(VR) 2 - ^(VS) 2 + l -AR + AVS - l -A 2 - . 

At the same time, the real- valued functions 9 = exp(R + S) and 9* = exp(R — S) 
obey the adjoint parabolic equations (to be compared with (35), (37)): 

8A = \A9* + - c S jP* (A3) 

d t 9 — — ~A0 + V(A0) + cj,0 



2 

where 
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c^ = -VA--A + 2Q-0 (A4) 



2 2 

,1/2 



Art 1 / 2 

2Q = = (Vi?) 2 4 A/ t > 



Let us furthermore assume that dAi/dxj and are locally Holder continuous 
and moreover that 

|A|<^i(x 2 + l) 1/2 (A5) 
\caJ <K 2 (x 2 + 1) 

for all i = 1,2, 3. Then, for a certain T which may be smaller that the previ- 
ously chosen T, there exists a unique fundamental solution k(y,s,x,t) defined for 
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all x,y G R 3 , but confined to a time interval < s < t < Tq. 

The solution has the following properties: 

(i) For any n G C C (R 3 ), 9^(x,t) = J R3 k(y, s,x,t)r](y)d 3 y is a regular solution of the 
forward equation (A. 3) such that /zm^^x, t) = 

(ii) k is a regular solutions of the forward equation (A. 3) in (x, t) variables and a 
regular solution of the backward equation (A. 3) in the (y, s) variables. 

(hi) J R 3 k(y, s,x,t)d 3 y < Ciexp(C 2 x 2 ) for some positive G\ and C 2 . 

(iv) k is strictly positive i.e. k > 0. 

(v) k obeys the Chapman-Kolmogorov equation j R3 k(y, s, z,u)k(z,u,x,t)d 3 z = 
k(y, s, x, t) for all s < u < t < T . 

In the above, properties (i) and the first part of (ii) are in fact a definition of 
a fundamental solution, P2| . The second property (ii) and the estimate (iii) are 
canonical properties of the fundamental solution, |62|. The property (iv) comes 
as a consequence of another estimate (Lemma 2 in Ref. |62| ): for every e 6 (0, t) and 
every i6i? 3 there exist constants A and /i such that k(y, s, x, t) > Aexp[— /i(a ? — y) 2 ] 
for all y G R 3 and all s G [0, t — e]. The Chapman-Kolmogorov equation (property 
(v)) follows from the uniqueness of positive solutions of the forward equation (A. 3), 

us 

Consequently, what needs to be proven is the uniqueness of the integral kernel k. 
To this end, let us assume that there are two fundamental solutions. Then, for any 
i] G C C (R 3 ) we have 



Ui(x,t)= hi(y,s,x,t)rj(y)d y (A6) 

for i = 1, 2, and Ui(x, t) are the regular solutions such that urx, 0) = r](x). 
Let us define u = ui — u\. Then, for k > C2 there holds: 

T r 

u(x,t)\exp(-k x 2 )d 3 xdt < (A.7) 



— I f f (ki + k 2 )(y, s, x,t)n(ylexp(—kQX 2 )d 3 yd 3 xdt . 

JO JR 3 JR? 

But we have J R3 ki(y, s, x, t)d 3 y < Ciexp(C2X 2 ) for i — 1,2. Consequently: 



/ / \u(x,t)\exp(—k x 2 )d 3 xdt<2\\r]\\ sup CiT ( exp[— (k — C 2 )x 2 ]d 3 x < 00 . 

^0 J R 3 J R3 

(A8) 

By Theorem AB2 of Ref. [S2|, in view of u(x, 0) = we here arrive at u(x,t) = 0. 
Since n was arbitrary, we conclude that ki = k 2 . The uniqueness has been proved. 
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One may argue about the seemingly spurious time interval restriction to To < T. 
This time domain limitation for the kernel k can be easily justified. Namely, let 
us consider the forward equation (cf. A. 3) d t u = hAu + AS/u + cu and define 

-*2 

v = uexp[— g(x, t)\ where g(x,t) = and a > 0. The newly introduced function 
v satisfies: 

d t v = -Av + (BW)v + c v (A.9) 

with B = A+ \Vg and c = . + ~(V ' g) 2 + \Ag + (AV)g - d t g. Because of 
d t g = ' by taking a large enough we can make the equation (A.9) "dissipative" , 

|67| , and so construct its fundamental solution ko- A fundamental solution of the 



forward equation (A. 3) is the given by: 

k(y, s, x, t) = k (y, s, x, t)exp[g(x, t) - g(y, s)] . (A10) 

However if a is large, then t can run only a small interval [0, 1/a]. 

All the above features of the fundamental solution, (i)-(v), can be explicitly 
checked (by inspection) for the kernel (46), associated with the parabolic pair (45) 
in a constant magnetic field. They hold true for times 0<s<t<To<7r, except 
for the property (iii) which is true only if T < | . Indeed: 

27rs in(t - s) Jtf exp {~ COt ^2 ~ Vl ^ 2 + ^ 2 _ V^) 2 }} ex P(- x 2Vi + x 1 y 2 )dy 1 dy 2 

(All) 

= cos(t — s) exp[ [x x + x 2 )\ ■ 

Since [2n(t — s)]^ 1 ^ 2 J R exp[—(x 3 — y 3 ) 2 /2(t — s)]dy 3 = 1, to achieve (iii) it suffices 
to put Ci = 1 and C 2 = ^tanT , provided T < |. 

The uniqueness of the kernel for time T < | is proved along the same lines as 
this for the fundamental solution of (A. 3). 
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